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Abstract 

Accuracy in complex dielectric permittivity calculations in binary dielectric mix- 
tures in two-dimensions are reported by taking into account the shape of the inclu- 
sion phase. The dielectric permittivity of the mixtures were calculated using the 
finite element, and the permittivities were estimated by two different procedures. 
The results were compared with those of analytical models based on mean field ap- 
proximation and regular arrangement of disks. We have approached the problem 
emphasizing the finite-size behavior in which regular polygons with n sides were 
assumed to mimic the disk inclusion phase. It was found that at low concentrations, 
< 30 %, considering an decagon (n = 10) cause an error of < 0.1 % in the effective 
medium quantities compared with the results obtained using the analytical models. 
Key Words: Dielectric mixtures, composite materials, the finite element method, 
finite size scaling. 

1. Introduction 

To predict and to better design (tailor) composite materials for electrical applications, 
such as composite insulators [1, 2] and elecromagetic shields [3], etc, have been a chal- 
lenge for both theoretical and practical importance [4] . In early days of electromagnetics 
theory, effective electrical properties, i.e., conductivity, a, and permittivity, e, of sys- 
tems composed of two phases have been calculated analytically by effective mean field 
approaches [5, 6, 7, 8, 9] (ema) and regular arrangement of disks [10, 11] (rad). And 
nowadays, computer simulations have become as an alternative way of doing science closer 
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to experiment than theory but complenientory to both. Moreover, with the help of new 
computation techniques to solve partial differential equations, numerical simulations of 
more complex systems, such as systems with several components with arbitrary shapes, 
can be considered and desired properties can be calculated. However, like experiments, 
computer simulations produce data rather than theories and should be judged on the 
quality of those data. Accodingly, the reliability of the applied technique and accuracy 
of the obtained results must be checked and be verified either by analytical solutions or 
by performing the same calculations with a different technique or tool [12]. The latter 
procedure can be, for example, applications of the finite element and the finite difference 
methods to the same problem. The former one, to verify with an analyical solution, is 
not, on the other hand, an easy task in most of the cases, in which either there are no 
analytical solutions or approximations are considered in the analytical solutions. More- 
over, when the numerical results are taken into consideration, there are plausible sources 
of errors which originate from the model. The numerical errors can be eliminated by 
changing the order of the applied model or the discritization method used. 

In this paper, the accuracy in the numerical calculations of effective electrical proper- 
ties of a binary dielectric mixture is reported by taking into account the finite-size scaling. 
A field simulation software, based on the finite element method (fem), has been used, 
and the effective properties has been calculated by two different ways. The results were 
compared with those of two anaytical solutions, which are based on ema and RAD. 

2. Electrical properties of binary mixtures 

Electrical properties of materials can be described by the dependence of either the complex 
dielectric permittivity, e(w), or complex (ac) conductivity, <;{u}), on the frequency, w and 
on the external variables such as temperature, pressure, humidity etc. Both of these 
quantities can be expressed in terms of one and other. Therefore, a general complex 
dielectric response of materials can be describe with the complex dielectric susceptibility, 

e{Lu) = e + x'{Lo)-ix"{Lo) + ^- (1) 

<^(w) = ieo^e{u)) (2) 

where i = \f—\ and eo is the dielectric permittivity of free space, I/SGtt nF/m. e and a 
are the dielectric constant at optical frequencies and ohmic conductivity of the material, 
respectively. Moreover, the simplest form of dielectric relaxation, x(w), is oberved for 
dilute solutions (mixtures) and ferroelectric materials [13], and is expressed in Debye 
form [14] 

XH = x(0)[1 + *c.t]-i (3) 

where x(0) ^-iid t are the dielectric strength and relaxation time (inverse relaxation rate) 
of the polarization. At frequencies much higher and lower than inverse relaxation time, 
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Figure 1. (a) ema and (b) RAD approaches to dielectric mixtures. 



T ^, there are only three material quantities that explicate electrical properties, e, x(0) 
and a, 

s{lu) = e + 






for oj 3> T 



e(c.) 



X(0) 



iequj for cj < T ^ 

When the analytical solutions of binary mixtures in two-dimensions are considered, 
the EMA approach supposes two concentric dielectric disks with dielectric permittivities 
ei and £2 such that they are embedded inside the effective medium with dielectric per- 
mittivity, e, as presented in Fig. la. In the two-dimensional RAD approach, on the other 
hand, the composite medium is assumed to be composed of monodispersed inclusions 
(disks) at square lattice sides, as illustrated in Fig. lb. The effective complex dielectric 
permittivities calculated by ema [5, 6, 9] and RAD [10, 11], then, yield 



£2 £1 + £1 + £2 £1 <? - £1 q 



£ema — - 


£2 + £1 - £2 9 + £l g 




TT -TTKq + AK^q{A + B) 


£rass 


^^TT + TTKq + AK^q{A + B) 


A = 


£1 - £2 



(4) 



£1 +£2 



(5) 
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where, q is the concentration of the inclusion phases, which are denoted by subscript (2), 
(0 < (J < !)• The parameters A and B are functions of the radius of the inclusion phase. 



A 



2r3^ 



1 



oo oo 

r4 - 16m4 ^ 2^ 2^ 

m— 1 n— 1 m— 1 

r + 2m 



r — 2m 



(r — 2m)^ — 4n2 
r — 2?7i + 1 



(r + 2m)2-4n2 (r - 2m + 1)2 - (2n - 1)2 
r + 2771 — 1 



{r + 2m- 1)2- (2n- 1)= 



(6) 



B = 2r3 V ^- 

^-^ r'^ ^ (2m - 1 

m— 1 

oo oo 



(2m- 1)4 
r — 2m, 



n—l ni—l 

r + 2m 



(r — 2m,)2 — (n — 1)2 

r — 2m + 1 



[r + 2m)2 - (n - 1)2 (r - 2m + 1)2 - An? 
r + 2m — 1 



(r + 2m - 1)2 - 4n2 
j4 and B values converge quickly a constant value for n = m> 10. 



(7) 



3. Numerical calculations 

Analytical calculations of electromagnetic problems using Maxwell's equations, are lim- 
ited to geometrical constraints. For some simple geometries with a small number of mate- 
rials (regions) and symmetries, analytical solutions can be found [15, 16, 11, 17, 18]. The 
analytical solutions are obtained using methods of images [19, 20], orthogonal functions 
(Green functions) [21] and complex variable techniques (conformal mapping) [22, 20, 15]. 
The conformal mapping can only be applied to two-dimensional problems in which the 
third spatial axis is neglected. For more complex geometries and non-homogeneous re- 
gions composed of several materials, numerical solutions of partial differential equations 
and of integral equations have been developed [23] e.g., the finite difference method, the 
finite element method, the method of moments and the boundary element method. 

Numerical solutions of electrostatic problems within a non-conducting medium are 
based on solving Poisson's equation 



V • (eeoV(/)) 



-P 



(8) 



where (j), and p denote the electrical potential and the total charge in the considered 
region, respectively. Moreover, if the medium is conductive where no free charges and 
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Figure 2. The unit cell of square lattice with corners at ABCD used in the calculations. Dark 
region is the disk inclusion with £2 and (72 and the lighter region is the matrix media with ei 
and (71. 
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Table 1. Composite electric properties obtained from the analytical formulae. 



Model 


q 


£m 


£h 


<T/e, [S/F] 


EMA 
RAD 


0.1 
0.1 


2.28571429 
2.28586956 


2.42552966 
2.42602610 


0.13749080 
0.13752067 


EMA 
RAD 


0.2 
0.2 


2.61538462 
2.61743641 


2.95231624 
2.95942150 


0.16802403 
0.16845390 


EMA 
RAD 


0.3 
0.3 


3.00000000 
3.01008959 


3.62134985 
3.66005320 


0.20703632 
0.20939659 



sources of charges are allowed, then, the solution is given by 

V • (crV0) == 



(9) 



When the medium is a mixture of these two cases (lossy dielectric), it consists of dielectric 
and conductive components. The solution, then, becomes time dependent and is given 
by a complex electric potential in the region with the coupling of Eqs. (8) and (9), which 
is also known as the continuity equation. 



V • (crV0) + V 



I ("0^^) 







or equivalently in Fourier-space with frequency dependent properties, 

\/ ■{[i€ae{Lj)uj]\/<j)} = 



(10) 



(11) 



where no free charges are allowed in the region, due to conductivity of the medium (lossy 
dielectric). Note that e{u!) in Eq. (11) is given in Eq. (1). 

In this work, we have used a field calculation software, 'Ace', [24] based on the fem. 
A square lattice unit-cell with a hard-disk inclusion was assumed, as shown in Fig. 2. 
The boundary conditions were chosen as follows; along line [AB] was the ground potential 
level, T^ = V, along line [CD] it was 1 V, and the lines [AD] and [BC] were the symmetry 
lines (axes of reflection). The calculations were performed under steady-state periodic 
conditions. The region of interest is meshed using an triangular meshing technique in 
which we have limited number of triangles by ^ 8000 elements. The minimum triangle 
size was selected using the size of the considered inclusion geometry in the meshing 
procedure. Moreover, the quadratic shape function was used to solve in the fem. 

The complex permittivity of a heterogeneous medium can be calculated in several 
ways, e.g., (i) by using the total current density, j, and the phase difference, 9 [25, 26, 27, 
28], (ii) Gauss' law and losses [26, 29, 30] and lastly (iii) by using the average values of 
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Figure 3. Comparison of tlie two analytical models. The solid ( ), dashed (- - -) and 

dotted ( ) lines represent the ratios of the high frequency dielectric permittivity, em, low 

frequency permittivity, e^x, and the conductivity, cr/eo, values, respectively. The region above 
the thin chain (- • -) line marks the 1 %. 



dielectric displacement (D) and electric field (E) [26, 31]. We have used the first two ways 
(i & ii) to calculate the complex dielectric permittivities, s{u!) of the structures considered. 
The phase parameters, e and a, were frequency and voltage independent. The values of e 
and a are chosen such that the interfacial polarization observed has a relaxation time, r, 
around 1 s. This is achieved when the matrix phase has ei = ei = 2 and cri — 1 pS/m, and 
the inclusion phase has 62 = £2 — 10 and a2 — 100 pS/m. The normalized conductivity 
values, c/eo, are cri/eo — 0.113 S/F and CT2/eo ~ 11-3 S/F, respectively We have focused 
on two frequencies, 27r|j.Hz and 27rMHz which are, respectively, used as subscripts ((o.) and 
(M) for the appropriate dielectric permittivity values. At these frequencies, the influence 
of interfacial polarization was negligible, | loga;| 3> t^^, and the composite medium can 
be expressed with three parameters, £m, £|j. [= £m + x(0)] and cr/eo- 

4. Results and discussion 



The concentration dependence of the electrical properties calculated using the analytical 
formulae of Eqs. (4) and (5) differ at some concentration level. This is ilustrated in Fig. 3 
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by the ratio of resulting effective parameters. For concentration values higher than 30 % 
{q > 0.3), the effective medium quantities obtained from the two analytical approaches 
start to differ. The behavior of ratios for the dielectric permittivity at low frequencies, 
£^ and the conductivity, ct, were similar and the change with respect to concentration, 
q, was steeper compared to the ratio of the dielectric permittivity at low frequencies, 
£m- The difference between the models is due to the approximations and simplifications 
considered in the geometries. In RAD assumptions there are neighboring inclusions whose 
charge distributions (polarization) influence the polarization of the individual inclusions, 
however, in ema the polarization of the inclusion is only due to the interface between 
the inclusion and the matrix phases. Accordingly, in the numerical simulations, three 
concentration levels were selected, q = {0.1,0.2,0.3}, and the £m-j £|x- and cr/eo-values 
were calculated. The same quantities obtained from ema and RAD are listed in Table 1 
for comparison. The differences in the Em and cr/eo values calculated from the analytical 
models are in the fourth, third and second number after the decimal point for q = 
0.1, 0.2, 0.3, respectively. However, the change in e^ is larger for both models at the same 
concentrations. 

In the simulations, we have assumed that the inclusion phase was a two-dimensional 
object, a re^M^ar polygon with n sides, as displayed in Fig. 4. The polygons were generated 
using a circle with radius, r, and a contraint on the area of the polygons, q. Then, the 
radius, r, as a function of n and q is expressed as, 



nsin(7r/n) cos(7r/n) 



(12) 



The denominator inside the square root approaches tt as n -^ oo. The size of radius 
was also used for the meshing procedure of the computation domain where r„/15 was 
the size of the minimum triangle. Furthermore, this approach leads to the finite-size 
scaling considerations [32, 33] in which as n — > oo, the inclusion phase is a perfect disk. 
As mentioned previously, two different methods were used to calculate the electrical 
quantities of the binary mixture, and Table 2 presents the results. The first remark was 
that the difference between the obtained e^-values using two different approaches. The 
discripancy between ef^ and e^ increased as the concentration level, q, was increased. 
Moreover, when the values were compared top those of Table 1, except e^, the other 
calculated quantities had good agreement with the analytical models. 

Obtained quantities, Sm-, £|a.- and cr/eo, can be described by a trivial relation, which 
considers the finite-size behavior, 

/(n, q)-\^ai{n- 2)"i(«) + a2{n - 2)"2(9) (13) 

where A — f{oo,q). We have also scaled the above equation with n — 2 since the cal- 
culations were performed in two-dimensions. The A-values are presented in Table 2 as 
n ^ oo, except e^ all others were close to the data in Table 1. In Fig. 5, an example of 
the finite-size behavior is shown. A critical number of sides, ric, is defined, such that over 



TUNCER 



Table 2. Electrical parameters calculated by the fem. 



n 






a lei [S/F] 


^b 

^M 


e° 


a/sj? [S/F] 


<7 = 0.1 


3 


2.31830984 


2.55431772 


0.14540117 


2.31831000 


2.57482200 


0.14462477 


4 


2.29814032 


2.46885006 


0.14011433 


2.29814046 


2.48120081 


0.13990578 


5 


2.29154397 


2.44497599 


0.13866220 


2.29154405 


2.45548560 


0.13852250 


6 


2.28903383 


2.43636313 


0.13814092 


2.28903405 


2.44625458 


0.13804715 


7 


2.28780142 


2.43222046 


0.13789067 


2.28780167 


2.44182314 


0.13782876 


8 


2.28712945 


2.43000429 


0.13775702 


2.28712965 


2.43945636 


0.13769620 


9 


2.28672438 


2.42867095 


0.13767662 


2.28672448 


2.43803276 


0.13762992 


10 


2.28646432 


2.42782513 


0.13762567 


2.28646448 


2.43713033 


0.13759796 


11 


2.28629434 


2.42727218 


0.13759236 


2.28629442 


2.43654060 


0.13756319 


12 


2.28617467 


2.42688644 


0.13756913 


2.28617471 


2.43612935 


0.13754747 


14 


2.28602947 


2.42641794 


0.13754093 


2.28602973 


2.43562974 


0.13752268 


16 


2.28594550 


2.42615013 


0.13752482 


2.28594567 


2.43534447 


0.13750994 


18 


2.28589649 


2.42599233 


0.13751533 


2.28589650 


2.43517637 


0.13750239 


99 


2.28577330 


2.42560592 


0.13749210 


2.28577345 


2.43476492 


0.13748560 


oo 

« — n 9 


2.28576594 


2.42558544 


0.13749087 


2.28576605 


2.43474329 


0.13747696 


3 


2.69416747 


3.30221747 


0.18974833 


2.69416760 


3.36013939 


0.18911785 


4 


2.64801022 


3.07703286 


0.17563592 


2.64801032 


3.11023106 


0.17541128 


5 


2.62892034 


3.00133543 


0.17099827 


2.62892033 


3.02810579 


0.17088012 


6 


2.62304978 


2.97945868 


0.16966607 


2.62304982 


3.00451441 


0.16958761 


7 


2.62025295 


2.96934543 


0.16905189 


2.62025298 


2.99363860 


0.16899605 


8 


2.61869852 


2.96379171 


0.16871499 


2.61869871 


2.98767247 


0.16866931 


9 


2.61776224 


2.96046863 


0.16851354 


2.61776238 


2.98410506 


0.16847596 


10 


2.61716591 


2.95837257 


0.16838656 


2.61716598 


2.98185662 


0.16836172 


11 


2.61677916 


2.95701332 


0.16830424 


2.61677923 


2.98039874 


0.16827771 


12 


2.61650596 


2.95606556 


0.16824686 


2.61650611 


2.97938287 


0.16822349 


14 


2.61616711 


2.95488967 


0.16817569 


2.61616735 


2.97812228 


0.16815940 


16 


2.61598326 


2.95424269 


0.16813652 


2.61598338 


2.97742871 


0.16812236 


99 


2.61560667 


2.95293358 


0.16805731 


2.61560681 


2.97602608 


0.16806089 


oo 
n — n "? 


2.61561521 


2.95297927 


0.16806011 


2.61561540 


2.97607565 


0.16805491 


y — u.o 

3 


3.15432863 


4.44992747 


0.25956453 


3.15432877 


4.59647254 


0.25812644 


4 


3.07199643 


3.93508366 


0.22640982 


3.07199661 


4.00935583 


0.22583378 


5 


3.02527912 


3.72286805 


0.21324951 


3.02527918 


3.77630766 


0.21302790 


6 


3.01422774 


3.67689028 


0.21042426 


3.01422781 


3.72627694 


0.21031765 


7 


3.00968689 


3.65880191 


0.20931774 


3.00968704 


3.70668244 


0.20923851 


8 


3.00700955 


3.64824288 


0.20867221 


3.00700972 


3.69525097 


0.20859681 


9 


3.00525073 


3.64134674 


0.20825088 


3.00525070 


3.68779013 


0.20820236 


10 


3.00420698 


3.63728448 


0.20800287 


3.00420720 


3.68339788 


0.20795785 


11 


3.00352151 


3.63462285 


0.20784042 


3.00352166 


3.68052130 


0.20781141 


12 


3.00303949 


3.63276735 


0.20772723 


3.00303968 


3.67851699 


0.20769857 


14 


3.00245464 


3.63050408 


0.20758916 


3.00245485 


3.67607184 


0.20756988 


16 


3.00211061 


3.62921448 


0.20751059 


3.00211061 


3.67468077 


0.20749301 


99 


3.00143716 


3.62665177 


0.20735438 


3.00143741 


3.67191430 


0.20734844 


oo 


3.00147224 


3.62677466 


0.20736184 


3.00147232 


3.67204654 


0.20735887 



Calculated from the current and phase shift between the applied voltage and current. 
^ Calculated using Gauss' law and the total losses in the medium. 
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Table 3. Finite-size behavior modeling parameters. 





9 


logaf 


a\ 


log 4 


"2 


log a? 


-\ 


logfla 


a\ 


Em 
cr/eo 


0.1 
0.1 
0.1 


-1.500 
-0.930 
-2.145 


-0.215 
-0.223 
-0.224 


-4.174 
-3.530 
-5.077 


-0.017 
-0.022 
-0.010 


-1.500 
-0.895 

-2.178 


-0.215 
-0.224 
-0.230 


-4.172 
-3.487 
-5.014 


-0.017 
-0.023 
-0.003 


o-/eo 


0.2 
0.2 
0.2 


-1.104 
-0.492 
-2.202 


-0.220 
-0.230 
-0.230 


-3.994 
-3.240 
-5.253 


-0.015 
-0.022 
-0.008 


-1.104 
-0.455 
-2.210 


-0.220 
-0.231 
-0.233 


-3.990 
-3.195 
-5.330 


-0.015 
-0.022 
-0.005 


£m 


0.3 
0.3 
0.3 


-0.806 
-0.127 
-1.829 


-0.227 
-0.241 
-0.242 


-3.724 
-2.970 
-4.951 


-0.015 
-0.021 
-0.011 


-0.806 
-0.083 
-1.840 


-0.227 
-0.242 
-0.243 


-3.722 
-2.926 
-4.993 


-0.015 
-0.021 
-0.009 



^ Obtained for values calculated by using the current and phase shift between the applied 

voltage and current. 
'' Obtained for values calculated by using Gauss' law and the total losses in the medium. 



this value, n > ric, the effective properties of medium with regular polygons as inclusions 
are approximately similar to those of a medium with disk shape inclusions. The anal- 
ysis showed that Uc is approximately 15 regardless the concentration levels considered, 
q < 0.3, and the error in the calculations is < 0.01 % for ric > 15, as displayed in Fig. 5. 
In fact, even an decagon {n — 10) can imitate a disk in which the error in the calculated 
electrical quantities is less than < 0.1 %. Furthermore, Eq. (13) is divided in two, 

In Table 3, the parameters of this behavior, Eq. (13), from a curve-fitting procedure are 
presented. All calculated quantities had similar behavior as in the Fig. 5. It was clear 
to us that the exponents were concentration independent, a ^ F{q). Moreover, although 
ai was constant for all considered quantities, a2 was dependent on quantities. Finally, 
there was no trivial relation between concentration and obtained a values. 

5. Conclusions 

Calculations of effective electrical properties of binary dielectric mixtures were used to 
evaluate number of sides of a regular polygon in order the considered polygon to imitate a 
disk in computer calculations. We assumed a meduim with inclusions as regular polygons 
with n sides, and calculated the electrical quantities, i.e., dielectric permittivity and ohmic 

10 
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conductivity by using the FEM, which were later compared to those of analytical formulae. 
In the simulations the concentration of the inclusion phase was constant and the shape of 
the inclusion was assumed to be regular polygons with n sides. The size of the polygon 
was used to control the fem discretization of the computational domain. The analytical 
models were based on ema and RAD, and two different procedures were used to estimate 
the effective properties. It was found that the procedure based on Gauss' law and the 
total losses in the medium was not that succesfuU as the other one based on the current 
and phase shift between the applied voltage and current. In order to find the polygon to 
mimic a disk, the finite-size scaling behavior was introduced by considering the number of 
sides, ri, in the regular polygons; in reality, a regular polygon with n sides becomes a disk 
as n ^ oo. It was found that for n > 15, there was no significant change in the effective 
electrical quantities of mixture for low concentrations, q < 0.3. Consequently, when 
the effective quantities of mixture with an regular decagon inclusion was compated with 
those obtained from the analytical formulae for a similar mixture with a disk inclusion, 
the percentage error between quantities yield less than 0.1 %. 
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Figure 4. Finite-size scaling of the inclusion sliapes (polygons). 
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Figure 5. Normalized dependence of incrimental iiigii frequency dielectric permittivity, [em — 
A|/A on number of regular polygon sides, n. The symbols open (Q) and filled (•) indicate the 
solutions obtained using the current density and phase shift between the applied voltage and 
current and using Gauss' law and the total losses in the medium, respectively. The solid lines 
( ) represents the fitted curves, ric is the critical side number for regular polygons. 
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